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I .  INTRODUCTION 


The  importance  of  the  shaped-charge  warhead  as  an  antitank  weapon  has 
justified  repeated  efforts  to  analyze  the  phenomena  associated  with  its 
performance.  The  conventional  shaped-charge  munition  consists  of  a  cylin¬ 
drical  charge  of  high  explosive  with  an  axisymmetr ic ,  conical  cavity  at  one 
end  into  which  is  inserted  a  conical,  metallic  shell.  On  detonation  of  the 
explosive  from  a  point  on  the  axis  of  symmetry  of  the  charge  at  the  end 
opposite  the  cavity,  a  pressure  wave  propagates  toward  the  conical  metal 
"liner,"  and  passes  over  it,  sequentially  collapsing  the  shell  and  generat¬ 
ing  a  high-velocity  stream  of  metal,  called  the  "jet"  and  a  lower  velocity, 
more  massive  body,  called  the  "slug."  The  jet  travels  along  the  axis  of 
symmetry  of  the  system,  in  the  direction  of  the  propagation  of  the  detona¬ 
tion  wave,  and  is  capable  of  penetrating  the  high-strength  steel  used  to 
protect  military  vehicles. 

Much  of  the  research  done  on  shaped-charge  physics  in  England  and  in 
the  United  States  during  the  Second  World  War  was  summarized  in  1948  in  an 
article  by  Garrett  Birkhoff,  Duncan  MacDougall,  Emerson  Pugh  and  Sir 

Geoffrey  Taylor. ^  These  authors  derived  equations  describing  the  jet 
generated  by  the  symmetrical  collapse  of  a  hollow,  ductile  wedge  under 
steady  flow  conditions.  Because  of  the  extremely-high  pressures  acting  on 
the  explosively-driven  metal  wedge,  the  wedge  material  was  treated  as  a 
hydrodynamically-perfect  fluid.  Although  not  strictly  applicable  to  the 
collapse  of  hollow  cones,  this  theory  was  used  to  predict  the  mass  and 
velocity  of  jets  from  conventional  shaped-charge  munitions  and  the  pre¬ 
dictions  were  in  good  agreement  with  the  experimental  data  for  the  first- 
formed  portions  of  these  jets. 

The  steady-state  theory  of  jet  formation  from  wedge-lined  charges  was 

2 

modified  in  1952  by  Emerson  Pugh,  Robert  Eichelberger ,  and  Norman  Rostoker 
to  account  for  the  variation  in  the  collapse  velocity  imparted  to  the  wall 
of  the  metal  shell  of  a  conventional  shaped  charge.  The  introduction  of 
variable-flow  input  conditions  into  the  steady  jet-formation  theory  resulted 
in  a  description  of  the  velocity  gradient  found  between  the  tip  and  tail  of 
a  conventional  shaped-charge  jet. 

Experimental  evidence  supporting  the  theory  of  Pugh,  Eichelberger  and 

3 

Rostoker  was  presented  m  1952  by  Eichelberger  and  Pugh.  By  recovering  the 
jet  and  "slug"  material  from  conventional  shaped  charges  with  conical,  steel 
liners  and  by  measuring  the  velocity  of  different  portions  of  the  jet,  the 
authors  were  able  to  show  that,  within  the  experimental  error,  the  jets 
produced  from  collapsing  cones  were  well-described  by  the  "unsteady"  jet 
formation  theory  valid  for  collapsing  wedges. 


L 


By  1954,  Eichelberger  had  refined  his  experimental  technique  suffi¬ 
ciently  to  detect  small  but  systematic  discrepancies  between  the  observed 
characteristics  of  the  jet  from  a  collapsing  cone  and  the  characteristics 
predicted  by  the  unsteady  jet  formation  theory  for  collapsing  wedges.  He 
showed  that  the  finite  time  required  to  accelerate  the  liner  to  its  final 
collapse  velocity,  neglected  in  the  usual  theoretical  treatment,  did  not 
explain  the  observed  discrepancies.  Another  possible  source  of  error  in  the 
theoretical  treatment,  neglect  of  the  convergence-velocity  gradient 


described  by  Sterne,"’  was  suggested,^  but  not  explored. 


With  the  advent  of  the  two-dimensional,  Lagrangian,  hydrodynamical 


computer  codes,  such  as  HEMP,  it  became  possible  to  compute  in  detail  the 
initial  motion  of  an  axisymmetr ic ,  shaped-charge  liner.  Unfortunately,  for 
most  devices,  the  material  distortion  occurring  in  the  jet  formation  region 
could  not  be  successfully  treated.  However,  these  calculations  emphasized 
the  progressive  thickening  of  the  shell  wall  which  occurs  during  the 
collapse  of  a  conical  shell.  This  "convergence  effect"  leads  to  a  velocity 
gradient  through  the  thickness  of  the  shell  which  influences  the  jet 
formation  process.  Such  convergence-velocity  gradients  were  described  by 
Sterne  in  1950  and  suggested  as  possible  perturbing  influences  on  the 
simple,  hydrodynamic,  jet-formation  theory  of  Birkhoff,  MacDougall,  Pugh, 

and  Taylor.^ 


In  this  report,  the  quantitative  effect  of  the  "convergence-velocity 
gradient"  on  the  characteristics  of  the  jets  formed  from  collapsing  cones  is 
evaluated  in  the  spirit  of  the  previous  analytical  jet-formation  models.  A 
repartition  of  the  mass  and  velocity  in  the  jet,  differing  from  the  predic¬ 
tions  of  the  classical  theory,  is  shown  to  occur  under  these  conditions. 


II.  THE  CONVERGENCE-VELOCITY  GRADIENT 


Suppose  that  a  fluid  shell  of  mass  m  per  unit  length,  with  external 
radius  r 2  and  internal  radius  r^,  executes  a  purely  radial  contraction.  If 

the  fluid  is  incompressible  and  mass  is  conserved,  it  follows  that 


-t—  (r2  -  r2)  =  0 

dt  2  l 


(2.1) 


and  the  quantity 


r2  -  r2 


(2.2) 


is  a  constant  of  the  motion.  Equation  (2.1)  implies  that 


r  r  =  r  r 
2  2  11 


(2.3) 


vSsS 


where  r  =  — *-  •  More  generally,  since  the  mass  of  the  shell  can  be 

1  dt 

subdivided  into  smaller  radial  regions,  Equation  (2.1)  implies 


•  =  r  r  , 

rr  2  2 


^  <  r  <  r2  . 


(2.4) 


If  r^/r^  as  1,  there  is  initially  a  nearly  linear  velocity  gradient  through 
the  thickness  of  the  shell,  with  the  highest  velocity  at  the  inner  surface. 

Sterne  characterizes  the  shell  motion  by  the  variable  x,  the  ratio  of 
r^  to  r^.  Then,  for  example, 


r.  -  r  V  x  -  1 

1  o 


(2.5) 


The  kinetic  energy  per  unit  length  of  the  cylinder  is  found  by  integrating 
the  kinetic  energy  density  per  unit  length  over  the  radial  extent  of  the 
cylinder.  The  incremental  kinetic  energy  per  unit  length,  dT,  due  to  the 
element  of  mass  dm,  located  between  r  and  r  +  dr  and  having  characteristic 
radial  velocity  r>  is 


1  •  2 

dT  =  -  r  dm. 


dm  =  2irordr 


(2.6) 


(2.7) 


where  o  is  the  mass  density  of  the  shell  material,  so 


purr “dr 


(2.8) 


Since  rr  is  constant  throughout  the  shell,  it  may  be  moved  outside  the 
integral  and 


T  =  tic  r  r 
2 


f2  p 

2  Jr.  r 


2  22  x“ 

r  x  x  In  — r - 

o  2  i 

x  -  1 


(2.9) 


(2.10) 


M  ’ 


This  kinetic  energy  is  assumed  to  remain  constant,  so  the  radial 
velocity  of  any  lagrangian  element  of  mass,  p,  may  be  found  from 


r  r 
2  2 


r  = 

u 


(2.11) 


where 


u  = 


to  (r2  -  r2  ) 

2  M 


(2.12) 


Then 


r  = 

V 


r  r 
2  2 


/  2  “ 
r  -  p/irp 
2 


(2.13) 


r  xx 
o 


J — 2 — 

'  Y  — 


x  -  p/m 


(2.14) 


But,  from  Equation  (2.10), 


f~2T 

V 


ro  In  [x2/(x2  -  1)] 


(2.15) 


so 


. . . . . . . (2.16) 

/in  {x^/(x 2  -  1)]  /  x2  -  V/m 

The  predictions  of  Sterne's  theory  for  the  radial  collapse  of  a 
cylindrical  shell  which  conserves  its  mechanical  energy  are  shown  in  Figure 
1.  The  copper  shell  has  an  initial  external  radius  of  25mm,  an  initial 
external  radial  velocity  of  1.5  km/s,  and  an  initial  internal  radius  of 
23mm.  The  solid  linesin  Figure  1  represent  the  velocities  of  the  inner,  I, 
middle,  M,  and  outer,  0,  parts  of  the  shell.  The  dotted  lines  indicate  the 
positions  of  these  Lagrangian  elements  as  the  collapse  parameter  x  changes. 
As  the  variable  x  approaches  unity,  the  radial  velocity  of  the  internal 
surface  approaches  infinity  and  the  velocity  of  all  other  elements  of  the 
shell  approaches  zero.  Sterne  points  out  that  the  rapid  concentration  of 
kinetic  energy  at  the  inner  portion  of  the  shell  requires  the  generation  of 
very  large  internal  pressures. 

Figure  2  shows  the  collapse-velocity  vectors  computed  by  the  HEMP  code 
for  various  Lagrangian  elements  of  an  explosively-driven  conical  shell  which 
is  collapsing  toward  its  axis  of  symmetry.  Only  one-half  of  the  cross- 
section  of  the  conical  segment  is  depicted.  The  values  of  the  pressure 
£ield  intensity  (in  GPa)  at  the  centers  of  the  HEMP  computational  cells  are 
also  given.  The  gradient  of  the  velocity  field  through  the  thickness  of  the 
shell  is  apparent  as  well  as  the  divergence  of  the  shell  boundaries  as  the 
flow  approaches  the  axis.  A  wedge  of  equivalent  initial  thickness  explo¬ 
sively-collapsed  shows  no  such  velocity  gradient.  The  covergence-velocity 
gradient  is  then  formed  when  the  flow  divergence  velocities  are  superimposed 
on  the  explosively  driven  collapse  velocity. 


3r  (u,t)_  [1 T 

at  'V  m 


RADI/ 


Figure  2.  Cross-secion  in  the  upper  half-plane  of  the  HEMP  computational  grid 
simulating  the  explosively-driven  collapse  of  a  conical  shell. 
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III.  THE  ASSUMPTIONS  AND  RESULTS  OF  THE  PREVIOUS  STEADY-STATE 

JET-FORMATION  THEORY 

The  original  analysis  of  jet  formation  from  a  collapsing  wedge,  the 
Birkhoff,  MacDougall,  Pugh,  and  Taylor  (BMPT)  theory,  was  based  on  the 
premise  that,  for  steady  flow  conditions,  the  process  was  equivalent  to  a 
diverted  flow  issuing  from  a  moving  collision  point.  This  process  is 
illustrated  in  Figure  3,  taken  from  the  BMPT  paper,  Reference  1.  Figure  3a 
shows  the  formation  of  a  jet  and  slug  from  a  wedge-shell  whose  sides  col¬ 
lapse  with  constant  velocity  V  ,  driven  by  the  detonation  of  a  charge  of 

high  explosive  that  was  in  contact  with  its  outer  surface.  The  solid  lines 
in  Figure  3a  show  the  collapse  conditions  at  one  instant  of  time  and  the 
dotted  lines  show  the  conditions  at  a  later  instant.  Figure  3b  shows  the 
formation  of  jet  and  slug  by  the  wedge-shell  shown  in  Figure  3a  from  the 
viewpoint  of  an  observer  stationed  at  the  moving  junction  A. 

The  high  pressures  involved  in  the  process  far  exceed  the  mechanical 
strength  of  the  metal  in  the  shell,  so  this  material  was  assumed  to  behave 
like  a  perfect  fluid.*  To  eliminate  the  complexity  of  pressure  waves  and 
because  the  metal  of  the  shell  is  not  very  compressible,  the  fluid  was 
considered  to  be  incompressible.  The  flow  of  material  through  the  collision 
region  was  determined  by  applying  the  laws  of  conservation  of  linear  momen¬ 
tum,  mechanical  energy,  and  mass  at  the  boundaries  of  the  collision  region. 

In  BMPT  theory,  the  constant  velocity  of  the  wedge  wall,  V  ,  is 

resolved  into  a  "flow"  component,  Vp,  along  the  wall  and  a  "stagnation 

point"  component,  V  ,  along  the  axis  of  symmetry  of  the  shell.  In  the 

sp 

reference  frame  of  the  collision  region,  traveling  down  the  axis  of  symmetry 
with  vsp,  the  collapse  velocity  is  the  same  along  each  streamline.**  In 

terms  of  the  collapse  angle,  /2  ,  and  the  y-  and  x-components  of  the  collapse 
velocity,  VQy  and  V q^,  respectively,  the  magnitude  of  the  flow  velocity  is 


V  =  V  esc  /3  , 
F  oy  r 


(3.1) 


while  the  stagnation  point  velocity  of  the  moving  junction  A  has  a  magnitude 
given  by 


V  =  V_  cos  B  +  V  . 

c  r\  h  '  r*  v 


(3.2) 


In  the  stationary  frame  of  reference,  the  jet  velocity  magnitude,  V  is 


V.  =  V_  +  V 
j  F  sp 


*A  perfect  fluid  is  non-viscous  and  non-conducting. 


(3.3) 


**A  streamline  maps  the  direction  of  the  velocity  at  every  point  in  space  at 
any  instant.  The  paths  of  the  fluid  particles  are  tangent  to  the  stream¬ 
lines  at  any  instant. 


since  the  jet  flow  "turns  the  corner"  and  proceeds  in  the  same  direction 
that  the  stagnation  point  is  moving.  The  magnitude  of  the  slug  velocity, 


V  -  V 
sp  F* 


(3.4) 


since  the  slug  flow  is  opposite  to  the  stagnation  point  motion.  The  mass 

partition  into  jet  mass,  m.,  and  slug  mass,  m  ,  of  the  incoming  liner 

J  3 


(shell)  mass,  m^,  is 


—  =  cos2  (B/2) 
m 

1 


(3.5) 


=  sin2  (B/2) 


(3.6) 


The  conclusions  reached  through  this  analysis  were  considered 
applicable  to  jets  produced  from  both  conical  shells  and  wedge  shells.  In 
the  absence  of  convergence  velocity  gradients,  the  same  results  can  be 
derived  for  conical  geometry,  by  distorting  the  geometry  somewhat  to  mask 
the  essentially  divergent  nature  of  the  flow.  In  the  next  section,  the 
equivalent  theory  for  a  conical  geometry  is  developed,  under  the  assumption 
of  "convergenceless  flow,"  in  order  to  lay  the  foundations  for  the  new 
theory  which  explicitly  addresses  the  convergence  effect. 


IV.  JETS  FROM  CONVERGENCELESS  CONE  COLLAPSE 

In  order  to  approximate  the  conditions  required  for  the  application  of 
the  steady-state,  BMPT,  wedge-jet  theory  to  a  collapsing  conical  shell,  the 
thickness  and  collapse  velocity  of  the  shell  are  assumed  to  be  adjusted  so 
that  the  collision  region  moves  with  constant  velocity  down  the  axis  of 
symmetry  of  the  cone  and  the  mass  entering  the  collision  region  is  constant. 
Then,  from  a  reference  frame  in  the  collision  region,  the  incoming  and  out¬ 
going  flows  are  steady.* 

The  conditions  of  the  liner  collapse  are  illustrated  schematically  in 
Figure  4.  Half  the  metal  liner  is  shown  in  cross-section  at  the  top  of  the 
figure  along  with  the  charge  of  high  explosive.  At  the  bottom  of  the 
figure,  the  conical  liner  is  shown  partially  collapsed  with  a  small  amount 
of  undetonated  explosive  still  remaining  at  the  end  of  the  charge. 


♦The  flow  is  steady  if  there  is  no  variation  with  time  in  either  the 
magnitude  or  the  direction  of  the  velocity  at  any  stationary  point  in  the 
space  through  which  the  fluid  flows. 
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Figure  4.  Schematic  diagram  of  "constant-velocity"  shaped  charge  with  tapered 
conical  shell  liner;  (a)  before  detonation  of  the  high  explosive, 

(b)  with  the  explosive  partially  detonated. 


The  collapse  conditions  are  such  that  the  collapse  angle  0 ,  measured 
with  respect  to  the  central  streamline  of  the  incoming  flow,  remains 
constant.  Since  an  individual  liner  segment  swells  as  it  approaches  the 
axis  of  symmetry,  the  incoming  flow  is  no  longer  uniform,*  as  it  was  in  the 
case  of  the  wedge  flow.  The  incoming  flow  is  divergent  as  it  approaches  P 
and  the  separating  streamlines  have  relative  velocities  with  respect  to  each 
other.  A  convergence  velocity  gradient  necessarily  exists  which  will  be 
neglected  in  the  analysis  of  this  section. 

Figure  5a  shows  the  areas  intercepted  by  the  flow  on  an  imaginary 
surface  enclosing  the  jet-formation  region.  Figure  5b  shows  the  boundaries 
of  a  "quasi-uniform"  flow  through  the  formation  region,  reminiscent  of  wedge 
flow.  The  incoming  flow  is  represented  as  parallel  flow  because  constant 
velocity  divergent  flow  raises  many  of  the  same  complications  associated 
with  the  actual  flow  to  be  treated  later.  The  adjustment  in  the  apparent 
"thickness"  of  the  flow  brought  about  by  the  convergence  of  the  mass  in  the 
conical  shell  and  the  consequent  rearrangement  of  the  spacing  between 
streamlines  occurs  in  this  "convergenceless"  model  in  the  interior  of  the 
collision  region  and  is  "invisible"  at  the  boundaries. 

In  Figure  5b,  let  CE  ■  a  be  the  initial  thickness  of  the  flow  as  it 
enters  through  the  bounding  surface.  Since  the  flow  is  axisymmetric ,  the 
total  area  of  the  boundary  surface  intersected  by  the  flow  is  the  area  of  a 
ribbon  of  width  a  and  length  2f  r, 


A  =  2irra  =  ir(r  _  +  r  )a  (4.1) 

o  u  t, 

on  the  same  diagram,  let  FG  *  a^ ,  and  BA  *  aj  be  the  radii  of  the  deflected 
flows.  The  boundary  surface  areas  intersected  by  these  flows  are 


Ai  =  it  a2  (4.2) 

1 

and 

a2  *  "a*  (4.3) 

In  the  reference  frame  of  the  moving  collision  region,  the  incoming  momentum 
per  unit  length  of  quasi-uniform  flow  has  a  z-component  given  by 


$a(E)  p  (v^cosB)  2 it  (r^-ex) 

o(c) 

-2iroVj,  cos  8  (rca  -  ea2/2)  , 


dx 


(4.4) 


♦The  flow  is  uniform  if  all  streamlines  remain  parallel  to  one  another. 
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where  x  is  a  coordinate  which  measures  the  thickness  of  the  flow  as 
indicated  in  Figure  5b, V  is  the  magnitude  of  the  flow  velocity. 


c  =  <rc  -  rE)/a, 


(4.5) 


p  is  the  mass  density  of  the  shell  material,  p  is  the  collapse  angle  defined 
by  an  incoming  streamline  and  the  axis  of  symmetry,  and  p^n  is  directed 

alonp  the  negative  z-axis. 

The  outgoing  momentum  per  unit  length  has  a  z-component  of  the  form: 


-r 

J  o 


a2(B) 


/•M(F) 

pV  2-rydy  +  |  pV  2try'dy' 

J  o(G)  F 


(4.6) 


where  y  and  y'  are  indicated  in  Figure  5b.  For  continuous,  incompressible 
flow,  the  relation  between  a,  a^  and  a^  is 

0”(r  +  r_)  a  =  o,na2+o'’Ta2 


c  'E 


2  1 


al  =  ^ (r  +  r  )  a  -  a2 
C  E  2 


P°U  =  -2ttPVf  [a2  -  (rc  *  r£)a/2]  . 


If  linear  momentum  is  conserved  in  the  collision  frame, 


cos  2  (r  a  -  ea2/2)  =  a2  -  (r  +  r  )  a/2 

C  2  C  E  2  ' 

a2  =  (r  +  r  )  ^  (cos  3+1) 

2  u  e  2 


(4.7) 


(4.8) 


(4.9) 


(r  +  r  )a  cos2  (2/2) 
C 


(4.10) 


The  fraction  of  the  mass  in  the  incoming  flow  which  goes  into  the 
’slug"  is  then 


m  a  (r  +  r  ) 

1  C  E 


=  cos2  (2/2)  , 


(4.11) 


which  is  the  same  result  obtained  in  BMPT  theory  for  the  slug  from  the 
collapse  of  a  wedge  shell.  The  mass  partition  in  the  jet  is 


sin2  (B/2). 


(4.12) 


The  velocity  of  the  slug  portion  of  the  mass  will  be  V  -  V_;  that  of  the 

Bp  F 

jet  mass  will  be  Vgp  ♦  V^,  as  indicated  in  Figure  3b. 

This  jet-formation  model,  which  ignores  convergence  effects,  gives  a 
good  qualitative  description  of  the  jets  produced  by  the  collapse  of  conical 
shells.  When  the  effect  of  a  time  varying  collapse  velocity  and  collapse 
angle  is  incorporated  by  simply  changing  the  input  conditions  to  the  steady- 
state  model,  the  description  is  so  good  that,  in  the  first  examination  of 
the  validity  of  the  model,  the  predicted  and  experimentally  -  determined  jet 
characteristics  were  shown  to  agree  within  the  precision  of  the  experimental 
3 

data  .  Later,  when  a  more  sensitive  analytic  approach  was  taken  using 

4 

refined  experimental  data  ,  small  systematic  deviations  between  predicted 
and  observed  properties  were  detected,  but  no  theoretical  explanation  for 
the  effects  was  offered.  In  the  next  section,  the  influence  of  a  conver¬ 
gence-type  velocity  gradient  on  the  mass  and  velocity  partition  in  wedge— jet 
formation  is  investigated.  Small,  but  systematic, deviations  from  the 
results  of  the  BMPT  theory  are  shown  to  occut.  The  derived  formulae, 
although  artificial,  illustrate  the  changes  in  the  jet  properties  with  the 
least  mathematical  complexity. 


V.  THE  INFLUENCE  OF  A  CONVERGENCE-TYPE  VELOCITY  GRADIENT  OB  THE  MASS 
AND  VELOCITY  PARTITION  IB  REFLECTED  WEDGE-FLOW 

The  effect  of  a  velocity  gradient  on  jet  formation  in  the  stagnation- 
point  frame  of  reference  is  most  easily  assessed  by  using  the  wedge-flow 
model  first  proposed  by  Birkhoff,  MacDougall,  Pugh  and  Taylor.  The 
advantage  of  this  model  is  its  particularly  simple  flow  pattern. 

In  Figure  6,  one  half  the  cross-section  through  the  flow  and  bounding 
surface  is  shown.  The  initial  thickness  of  the  wedge  flow  is  CE  ~  a,  and 
the  thicknesses  of  the  deflected  forward  and  rearward  flows  are  FG  -  a^  and 

BA  -  aj ,  respectively.  The  velocity  gradient  across  the  incoming  flow  iB 

assumed  to  be  linear,  in  particular 


v_  (x)  =  V  •+  ctx 


0<  x  <  a 


(5.1) 


where  V  is  the  magnitude  of  the  flow  velocity  at  position  x,  measured 
F 

across  the  thickness,  with  x  *  0  being  point  C  on  Figure  6.  V  is  the  flow 

OF 

velocity  at  C  and 


a  =  (V  -  V  )/a 

E  c 

is  the  "strength"  of  the  velocity  gradient. 


(5.2) 


Since  the  flow  is  assumed  to  be  uniform  and  steady,  there  is  a  smooth 
division  of  the  streamlines  in  the  collision  region  and  the  incoming 
velocity  gradient  is  preserved  in  the  exiting  flows. 

The  incoming  momentum,  through  a  surface  CE  across  the  cross-section  of 
the  flow  and  extending  for  unit  length  perpendicular  to  the  plane  of  Figure 
6,  has  an  average  z-component  given  by 


_  in 

<PZ  >  =  -oa  <vp>  cos  B, 


(5.3) 


where  <V  >  is  the  average  flow  speed*  defined  by 
r 

L 


<V„>  H 


1  A  VFN  dA 


(5.4) 


VpN  being  the  component  of  the  flow  velocity  normal  to  the  area  element  dA. 


Then,  since  the  surface  CE  is  normal  to  the  flow 


.  m 

<?2  >  =  -oa  cosB 


Pi 

J  o  a 


(x) 


dx 


=  -0  cosB  (V  a  +  aa2/2) 


(5.5) 


The  incoming  mass  per  unit  time  is  pa  <Vp>.  The  outgoing  average 
momentum  has  the  z-component, 


out 

<p  >  = 

rz 


-oa2  <V  >  +  Daj  <v  > 
F2  FI  ' 


(5.6) 


♦This  definition  is  given,  for  example,  by  Jerzy  A.  Owczarek  in  Introduction 
to  Fluid  Mechanics,  (International  Text  Book  Company,  Scranton,  1968),  p. 
llTT 


and 


Then  since  , 


fa2  V*'>  a 

<V  >  =  I  -  dx‘  =  V  +  — 

F2  Jo  a2  OF  2  2‘ 


F  (x)  „  ,  a 

-  dx  =  V  +  -r-  a, 

a  OF  2 


V  ai  +  aaia2  +  -^-a2=V  a  +  -2-  a2  -  v  a 
OF  2  1  OF  2  of  2 


(5.12) 


(5.13) 


(5.14) 


Since  momentum  is  conserved, 


a  cos  3  (V  +  aa/2)  =  2a  (V  +  aa  /2 )  -  a  (v  +  aa/2) 
OF  2  O-  ?  OF 


(5.15) 


The  thickness  of  the  outgoing  flow  is  then 


V2  +  aa  (cos  3+1)  (v  +  aa/2) 
Or  nr 


(5.16) 


where  the  limitation 


0  <  a  <  a 
-  2  - 


determines  the  root  to  be  chosen.  Note  that,  if  Equation  (5.16)  is  expanded 
by  means  of  the  binomial  theorem,  the  result  is 


a  =  a  cos2(S/2)  +  ~ —  cos"  (S/2)  +  ... 
2  OF 


(5.17) 


.  %  4V  v 

>  *V‘\- 


V. 


V. 


in  agreement  with  BMPT  theory  if  a  =  0.  For  constant  influx  energy  (see 

Appendix  A),  as  a  increases,  V  decreases  until  finally,  the  limiting  value 

OF 

for  VQF  non-negative  is  =  a  cos  (S/2) 


As  one  would  expect  from  a  velocity  gradient  of  the  sort  assumed  (with 
a  >  0 ) ,  less  mass  is  required  in  the  jet  to  balance  the  axial  momentum  in  the 
stagnation  point  frame  of  reference  because  the  average  flow  velocity  in  the 
jet,  <Vp^>,  is  greater  then  the  average  flow  velocity  in  the  slug,  <7^^* 

Some  examples  of  the  phenomenon  are  shown  in  Table  I. 


Table  I,  Dependence  of  Wedge-Flow  on  Velocity  Gradient  Tor  Constant 
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VI.  JETS  FROM  COLLAPSING  CONICAL  SHELLS 


The  steady  flow  associated  with  a  collapsing  conical  shell  is  non- 
uniform.  This  introduces  additional  complication  into  its  analysis. 

According  to  Sterne's  theory  of  the  collapsing  cylindrical  shell,  the 
radial  velocity  of  the  central  streamline  is  essentially  unaffected  by  the 
internal  pressures  until  the  collapse  is  nearly  complete.  Consequently,  the 
other  streamlines  diverge  from  this  streamline  and  the  central  streamline  is 
a  natural  point  of  reference  in  the  flow  analysis. 

In  the  following  it  is  assumed  that,  in  a  reference  frame  in  which  this 
streamline  is  at  rest,  the  flow  pattern  appears  stationary.  This  conical 
flow  pattern  is  illustrated  in  Figure  7  along  with  the  flow  velocities  as 
seen  in  the  "stagnation  point"  reference  frame  of  the  central  streamline. 

A  boundaTy  region,  indicated  by  the  dashed  line  in  Figure  7,  is 
established  just  outside  the  high  pressure  collision  region  whose  existence 
was  indicated  in  Figure  2.  The  flow  velocities  are  obtained  by  resolving 
the  collapse  velocities  along  their  respective  streamlines.  The  reference 
frame  velocity  is  the  axial  component  of  the  collapse  velocity  on  the 
central  streamline,  evaluated  at  the  point  where  the  streamline  first 
intercepts  the  surface  of  integration,  just  above  the  collision  region.  The 
collapse  velocities  of  the  other  streamlines  have  different  axial  compo¬ 
nents,  since  they  are  diverging  from  the  central  one. 

The  momentum  of  the  fluid  entering  the  dotted  boundary  surface, 
illustrated  in  Figure  7,  will  be  conserved.  As  in  the  previous  analysis 
of  conical  flow,  CE  *  a  is  the  thickness  of  the  incoming  flow  and  the 

projected  area  of  the  flow  on  the  boundary  surface  is  A  =  ir(r  +  r  )  a, 

O  0  E 

where  rc  and  r£  are  radii  to  the  boundary  points  C  and  E  in  Figure  7. 

The  radii  of  the  outgoing  "jet"  and  "slug"  flows  are  a^  and  a^ ,  respec¬ 
tively.  The  projected  areas  of  the  outgoing  flows  on  the  boundary 

surface  are  A  =  ira2,  for  the  jet,  and  A  =  ua2  for  the  slug. 

11  2  2 

Let  x  label  the  position  of  a  streamline  in  the  incoming  flow  and 

let  V^(x)  be  the  velocity  component  of  the  incoming  flow  at  x  normal 
to  the  boundary  surface.  It  is  assumed  that  the  velocity  gradient  of 
this  normal  flow  is  linear,  i.e., 


%'*> 


VA_,  +  ax 
OFN 


0  <  x  <  a 


(6.1) 


where  is  the  normal  component  of  the  flow  velocity  at  point  C 

in  Figure  7,  and  x  is  measured  from  point  C.  The  velocity  components 
parallel  to  the  boundary  CE  are  associated  with  the  "natural"  stream¬ 
line  divergence  and  are  neglected  in  the  analysis  of  the  collision. 


i 

i 

( 


The  average  momentum  per  unit  length  of  the  flow  entering  through 


the  boundary  surface  has  a  z-component ,  <n1T1>  given  by 

*  ■ *9  9 


<P^>  =  -o  <VFN>  Ac  cos  3, 


where 


<v  > 

FN 


-u: 


V_, (x)  2ir  (r  -  ex)  dx. 
FN  ~ 


(6.2) 


(6.3) 


The  integrated  incoming  momentum,  in  the  negative  z-direction,  is 


<pxn>  =  -2-rro  F(a)  cos  B  , 
z 


where  F(a)  =  rca  Vqf;,  +  (nrc  -  eV^)  a^/2  -ciEa3/3, 


(6.4) 

(6.5) 


The  average  momentum  leaving  the  bounding  surface  has  a  z-component, 

<p°ut>  »  given  by 
z 


t  fa  (B)  /-a  (F) 

<p  >  -  -  I  2  oV  (x ' ) 2itx ' dx '  +  I  1  oV  (x")  2Trx"dx"  ,  (6.6) 

2  -'o(A)  F  Jo(G)  F 

where  Vp(x'  )  and  VF(x"  )  are  the  streamline  flow  velocities  in  the  slug 

(x' )  and  jet  (x"  )  respectively.  These  flows  are  normal  to  the  boundary 
surfaces. 

It  is  assumed  that  the  normal  components  of  the  individual  streamlines 
are  preserved  in  this  "perfect"  collision.  The  spacing  between  the  stream¬ 
lines  changes  due  to  convergence  effects.  For  continuous  incompressible 
flow,  conservation  of  mass  demands  that 


where 


and 


+  V  V 

=  a2<V  >  -t-  a2<V  > 
1  F1  2 

2F  (a) 

%  v ,,,, ^  — 
FN 
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<V  >  =  —f 
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■  j  1  V _(x")x"dx" 
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ra2 

<V  >  =  ~ T 
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j  VF  (x’lx’dx' 

^  0 

(6.7) 

(6.8) 

(6.9) 


(6.10) 
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Then,  the  expression  for  continuity  is 


rai  ra 2 

F (a) =  I  Vp (x") x"dx"+  I  V 

Jo  Jo 

=  f  1  (x") x"dx' 


V  (x')x'dx' 
F 


=  F (a)  -  I  VF(x')x,dx’ , 


<V  >  =  - — 

FI  2 

a 
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The  axial  component  of  the  average  momentum  is  therefore 


out 

<D  > 
*  Z 


-2-c  S2  f 

I  J  c 


VF(x')x'dx'  -  F(a)V  . 


Since  momentum  is  conserved, 


'  (a) cos? 


VF(x'5x'dx'  -  F(a)  , 


j  (1  +  cosS)  F(a)  =  I  VF (x ' ) x 'dx ' 


I  =  I  V  (x')x'dx' , 

2  I  0  F 


<v  „> 

F2 


.-:X-XvV.  -y/.'WS: 


In  order  to  evaluate  the  integral  I ^ »  the  velocity  distribution  in  the 
exiting  flow,  Vp(x,)>  must  be  determined  from  the  spacing  between  stream¬ 
lines  in  the  slug  flow.  Figure  8  shows  the  flow  variables  considered. 


Let  a^  be  the  lagrangian  label  of  the  streamline  which  divides  the 
jet-slug  flow.  The  area  of  a  streamtube*  in  the  incoming  slug  flow  is  then 


Ay  =  ”y  (ra*  +  ry) ,  0  <  y  <  a*  (6.18: 

2  ~  -  2 

where  Av  is  the  area  between  the  dividing  streamline,  a*,  and  the  stream¬ 
line  with  coordinate  y;  r  *  is  the  radial  position  of  the  point  a*  on 

a2  ^ 
the  boundary  surface  and  ryis  the  radial  position  of  the  point  y  on  this 

same  surface.  Let  the  bounding  surface  CE  be  normal  to  the  central 

streamline  of  the  flow. 

If  d  is  the  angle  between  the  line  CE  and  the  positive  z-axis  in 
F igure  8  ,  then 


and  also 


so 


r  *  =  r  -  ra* 

a  C  o 

2 


(6.19) 


(6.20) 


(6.21  ) 


*A  streamtube  is  a  tube  in  the  flow  velocity  field  filled  with  the  flowing 
fluid  and  formed  by  all  the  streamlines  passing  through  a  closed  curve 
defining  the  circumference  of  a  cross-section  of  the  tube. 


r»  * 


Figure  8.  Flow  configuration  for  steady,  non-uniform  flow  for  evaluatinq 
streamline  spacing. 
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where  £  is  given  by  Equation  (4.5).  The  area  between  any  two  streamlines  is 
then 


A v  =  "y  (2r  -  2ea*  +  ey)  . 

y  C  2 


(6.22) 


All  these  streamlines  are  carried  into  the  slug,  where  the  spacing  is 

determined  bv  the  area  A  /  of  the  stream  tube,  where 

x 


.  .  =  "X 


Since  the  flow  is  incompressible,  A  ,  *  A  and,  since  y  is  positive, 

x  y 

"(rC  -  ea2)  Jirc  '  a2)2  +  Ex’2 

V  ~  ■  + - - - - - 


(6.23) 


(6.24) 


Mapping  the  streamline  velocities  from  one  configuration  of  the  flow, 

V  (v)  incoming,  to  the  other,  V  (x')  outgoing,  leads  to  the  trans format  ion 
r  '  r 


Letting  A(a*) 


Vv)  =  VF(a*>  -  ay  =  V0F,;  +  aa* 


r  -  ra*  it  follows  that 
c  i 


(6.251 


Vx>)  =  vofn  +  ^a* +  7  A(a;! 


2j  A2  (a*)  +  £x'" 


(6.26) 


& 


« 

t 

s 


Vi 
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Substituting  for  V  (x  )  in  Eq.  (6.17)  for  the  integral  I  and  integrating 
•  F  2 


give  8 


V  or 

T  OFN  o 

1  =  - T~  a  + 

2  2  2  2 


2  C  n  if, 

+  a' - -  ' 

2  2e  2  3e-  L 


(r  _  -  ea*)  2  +  ea1 

c  2  : 


1  7/2  +  _i!_ 

J  I,2 


(6.27) 


(r  -  ea*)  3 
c  2 


The  result  of  equating  the  incoming  and  outgoing  "slug"  areas  is  that 


r  -  Wr2  -  ea2 

C _ y  c _  -> 


1  2  1  9  ,  1  0  -2  a  1  ,  9  o .  -3  /o 

V  2  VQFNa2  +  2  "rca!/C  "  IarCA  +  3^T  l*C  ■  Ca2]  V 


(6.28) 


(6.29) 


Now  1^  is  just  the  difference  between  F(a),  Eq .  (6.5),  and  1^.  Once  1^ 

is  known,  <Vpj>  follows  from  Eq .  (6.13).  Then  <^VF2'>  can  be  obtained  bv 

dividing  twice  I0  by  a- .  The  magnitude  of  the  average  flow  velocitv  into 
the  slug  is  then 


-7s-  ~  <4  -  C4  -  **n  3/2>  • 


<V  >  =  V  + - - (ri  -  \rt  -  ea^ 

F2  CFN  e  3e2a2  0  L  c  2 

2 


(6.30) 


The  expression  for  conservation  of  the  axial  component  of  the  linear 
momentum,  Eq .  (6.16),  becomes 


1  a  2  2  3/2  arr 

(1  +  cos  2)  F(a)  =  —  (v  +  — a2  +  -  (r  -  ca  )  -  - 


2  OFM 


2  3  c 


.2  C  2' 


— -  (6.31) 

r  ^ 


j  (1  +  cosB)  F(a)  + 


M  =  2(VOFN  +  • 


(6.32) 


(6.33) 


(6.34) 


-  x*  ■J'V**_*  KT  m"  »  V  ‘-A  A  Jr  .% 


Then,  Equation  (6.31)  becomes 


K  -  Ma2  =  L(r2  -  ea?) 3/2 . 
2  C  2 


Squaring  both  sides,  rearranging  t«nns  and  letting 


u  =  a2 


M2  -  3<Lerc)2 

J7I2 


3L2er*  -  2 KM 


(6.35 


(6.36 


(6.37 


(6.38 


*’  -  Ac 


the  result  is  the  cubic  equation. 


■u3  +  Nu2  4  Vu  +  ?  *  0 


(6.39 


(6.40 


The  cubic  equation  for  the  mass  distribution,  Eq.  (6.40),  is  solved  by 
the  standard  methods.*  If 


3V—  N~ 


~  27p  -  2N3, 


■  {RW°3  *  r2) 


(6.41 


(6.42 


(6.43 


♦For  example,  see  Murray  R.  Spiegel,  Mathematical  Handbook  of  Formulas  and 
Tables ,  Schaum’s  outline  srries  (McGraw-Hill,  New  York,  1968)  ,  p. 32 . 


and 


then  the  solutions  to  Equation  (6.40)  are 


u  .  =  S  +  T  -  N/3  , 


u;  =  -  J  (S  +  T)  -  j  +  All  (S  -  T) 


(6.44) 


(6.45) 


(6.46) 


and 


u 

■5 


1 

2 


(S  +  T) 


N 

3 


(S  -  T). 


(6.47) 


i 


For  physical  reasons,-  u  a(r„  +  r^),and  this  criteria  is  used  to  select 
the  proper  root . 

Once  a,,  has  been  determined,  by  taking  the  positive  square  root  of  the 
appropriate  value  of  u,  ,  c^e  avera8e  fl°w  velocity  in  the  "slug"  flow, 

can  be  obtained  from  Eq .  ( 6 . 30 )  .<V^>,  the  average  "jet"  flow  velocity, 
follows  from  Eq .  (6.7),  since  a^  can  be  determined  from  the  area  identity, 


3(rC  +  V 


1  - 


a(rc  +  V 


(6.48) 


Table  2  shows  the  results  of  increasing  the  strength  of  the  velocity 
gradient,  n  ,  on  the  other  properties  of  the  flow  in  the  conical  geometry. 
For  these  conditions,  the  behavior  of  the  flow,  parameters  for  the  flow  with 
linear  convergence  -  velocity  gradient  are  quite  similar  in  character  to 
those  obtained  in  Section  5  for  wedge  flow  with  velocity  gradient,  The 
calculations  were  made  for  constant  influx  energy  (see  Appendix  C) .  For 
i  =  0,  the  treatment  of  Section  IV  is  used. 
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Table  2.  Dependence  of  Conical-Flow  Properties  in  Stagnation  Point 
Frame  on  Velocity  Gradient  for  Constant  Influx  Energy 


V0FN  -  2.900 

km/s 

a  ”  4.624 

linn 

,8-  36.4 

degree: 

£  -  0.824 

rc“  13'22 

mm 

rE  -  9.41 

mm 

a 

V 

0FN 

al 

a2 

<VFl> 

<V„„> 

F2 

-1  -1. 
i  mm  ) 

(kms  1 ) 

(mm) 

(mm) 

(kms  3 ) 

( kms 

0.0 

2.900 

3.19 

9.72 

2.900 

2.900 

0.1 

2.679 

3.06 

9.76 

3.116 

2.876 

0.2 

2.451 

2.97 

9.79 

3.323 

2.847 

0.3 

2.218 

2.86 

9.82 

3.548 

2.816 

0.4 

1.978 

2. B0 

9.84 

3.730 

2.779 

0.5 

1.732 

2.69 

9.87 

3.934 

2.740 

0.8 

1.479 

2.61 

9.89 

4  .149 

2.694 

0.8 

0.952 

2.46 

9.93 

4.520 

2.586 

1.0 

0.396 

2.33 

9.96 

4.881 

2.452 

VII.  COMPARISON  OF  THE  WEDGE  AND  CONICAL  SOLUTIONS 

Figure  9  illustrates  the  effect  of  varying  the  strength  of  the 
convergence-velocity  gradient,  a,  on  the  mass  partitioning  in  the  wedge  and 
conical  flows.  Both  flows  are  constrained  to  have  constant  influx 
thickness,  collapse  angle,  and  kinetic  energy.  For  the  wedge  flow,  the 
ratio  of  the  mass  of  the  slug  formed  by  a  collapsing  liner  element  to  the 

2 

mass  of  the  element  itself  increases  smoothly  from  cos  (/?/ 2)  to  cos  (/3/2). 
Above  cos  (/9/2),  because  of  the  constraint  of  constant  influx  energy,  the 
strength  of  the  gradient  can  only  be  increased  by  reversing  the  flow 
direction  of  part  of  the  flow,  a  physically  unappealing  situation.  For  the 
conical  flow,  the  fraction  of  the  mass  going  into  the  slug  increases 
smoothly,  following  the  wedge  flow  . 


Behavior  of  mass  partition  as  the  gradient  strength  increases  for  wedge  ai 
flow  under  the  condition  of  constant  influx  energy.  (Triangles  represent 


The  physical  reason  for  this  great  similarity  is  illustrated  with  the 
help  of  Figure  10.  The  quantity  F(x)  is 

r  V  x  +  (rsr  -  eV  )  x2/2  -  sax3/ 3 
C  OFN  C  OFN' 

and  plays  the  role  of  an  area-weighted  average  of  the  incoming  velocity.  As 
the  thickness  coordinate,  x,  is  increased  from  0  to  a,  4.624  mm,  the  average 
velocity  increases.  As  the  strength  of  the  velocity  gradient  increases,  the 
shape  of  the  velocity  distribution  changes  from  a  curve  with  gradually 
decreasing  slope,  to  one  of  constant  slope,  to  one  of  increasing  slope.  Even 
for  the  larger  values  of  a,  under  the  constraints  of  constant  influx  energy  and 

positive  V._It  the  contribution  from  the  cubic  term  is  not  great  and  the  distri- 
OrN 

bution  resembles  the  corresponding  wedge  distribution  quantity,  f(x)  =  V  x  + 
ax2/2. 

VIII.  APPLICATION  OF  THE  THEORY 

In  order  to  illustrate  the  new  jet-formation  analysis,  the  input 
boundary  conditions  from  a  HEMP  code  simulation  of  cone  collapse  will  be 
used  to  calculate  the  properties  of  the  jet  and  slug  segments  produced  from 
one  element  of  the  liner.  In  Figure  2,  the  material  in  the  segment  of  the 
cone  whose  cross-section  in  the  upper  half-plane  is  bounded  by  the  grid 
lines  8,  9,  14,  and  16,  is  assumed  to  collapse  steadily  into  the 
jet-formation  region  below  boundary  line  9.  The  velocities  and  positions  of 
the  node  points  on  the  boundary  surface  just  above  the  high-pressure 
collision  region  are  given  in  Table  3. 

Table  4  gives  the  flow  velocities  and  directions,  relative  to  the 
positive  z-axis,  of  these  boundary  points  as  well  as  the  "stagnation  point" 
velocities  for  all  three  streamlines,  as  computed  from  Equations  (3.1)  and 
(3.2),  with  the  appropriate  changes  in  notation  to  reflect  the  new 
coordinate  axes. 

The  "thickness"  of  the  incoming  flow,  as  measured  normal  to  the  central 
streamline,  is 

a  “  4.624  mm 

The  slope  of  the  boundary  line  in  the  flow,  as  defined  by  Equation  (4.5),  is 

E -  0.8240, 

where  the  radii  to  the  edges  of  the  boundary  surface  are 

rc  -  R(9 , 14)  -  13.22  mm 

and 


Tg  *  R(9,16)  ■  9.41  mm. 
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Table  3.  Input  Boundary  Conditions  For  The  Jet-Formation  Theory 


Node  V 

V 

V 

R 

Z 

z 

r 

Point  (km/s) 

( km/ s ) 

( km/s ) 

(mm) 

(mm) 

(9,14)  0.915 

-1.478 

1.738 

13.22 

-28.22 

(9,15)  1.153 

-1.803 

2.140 

11.40 

-27.06 

(9,16)  1.565 

-2.185 

2.688 

9.41 

-25.60 

4.  Flow  and  Stagnation  Point 

Velocities 

Computed  From  The 

Conditions 

Node 

V_ 

9V 

V 

F 

F 

sp 

Point 

(km/s) 

( deg) 

( km/ s ) 

(9,14) 

2.758 

32.4 

3.244 

(9,15) 

3.038 

36.4 

3.598 

(9,16) 

3.278 

41.8 

4.009 

In  order  to  analyze  this  data  with  the  BMPT  theory,  the  velocity 
gradient  and  flow  divergence  must  be  ignored.  Taking  the  central  streamline 
to  characterize  the  flow,  the  fraction  of  the  liner  mass  going  into  the 
"slug"  is,  from  Equation  (4.11), 


a(rc  +  rs) 


=  COS'  (2/2)  =  0 . 90?4 , 


since  P  =  36.4  degrees.  The  radius  of  the  slug,  a^ ,  is  9.72  mm.  The  jet 
velocity,  from  Equation  (3.3),  is  3.038  km/s  +  3.598  km/s  or  6.636  km/s, 
while  the  slug  velocity,  from  Equation  (3.4),  is  0.560  km/s. 

Note  that,  in  this  case,  the  incoming  kinetic  energy,  in  the  stagnation 
point  reference  frame,  is  expected  to  be  135.4  kJ/cm  (if  the  flow  had  no 
velocity  gradient). 


^  «*•  »*•  »  *  > 
■V  -.'.V.v  .V 


To  analyze  Che  data  with  the  new  theory,  the  (linear)  velocity  gradient 
across  the  thickness  a  must  be  determined.  The  components  of  the  flow 
velocities  of  the  surface  streamlines  normal  to  the  boundary  surface  are 

ICn  (9,14)1  =  |v*  (9,14)  [cos  4°  =  2.751  ^ 

and  _ _  . 

[  VpN  (9, 16)  I  =  |vF  (9 , 1*5)  1  cos  5.4°=  3.263  ^  . 


The  components  of  the  flow  velocities  parallel  to  the  boundary  surface  are 

[71^(9,14)1  =  -  !  V* (9,14)1  sin  4°  =  -  0.192  ~ 

i  FI  I  i  F  1  S 

and  , 

I  (9,16)  j  =  |  Vp  (9, 16)  [  sin  5.4°  =  +  0.308  — 


These  oppositely-directed  velocity  components  are  a  manifestation  of  the 
flow  divergence.  These  "lateral"  velocities  are  ignored  in  the  jet 
formation  calculation. 

The  strength  of  the  (normal)  velocity  gradient  is,  from  Equation  (5.2), 

a  =  0.1107 

mm 

Based  on  the  assumption  of  a  linear  velocity  gradient,  the  predicted  flow 
velocity  on  the  central  streamline  is  2.990  km/s,  about  98  percent  of  the 
input  value  of  3.038  km/s. 

The  fact  that  the  "stagnation  point"  velocities  of  the  surface 
streamlines  differ  from  each  other  and  from  the  stagnation  point  velocity  of 
the  central  (reference)  streamline  is  a  consequence  of  the  divergent  flow 
imposed  by  the  conical  geometry.  Nevertheless,  for  a  properly-tapered 
conical  shell,  these  streamlines  are  embedded  in  their  same  relative 
positions  in  the  surface  of  integration  as  this  surface  translates  down  the 
axis  of  symmetry  of  the  shell,  so  that  steady  flow  conditions  are  achieved. 

The  key  variables  in  the  theory  of  conical  flow  are  given  in  Table  5. 
The  cubic  equation  to  be  solved  is 

u’  +  Nu2  +  Vu  +  P  =  0 


From  Equations  (6.41)  and  (6.42), 


Q  =  -87.39  x  10  '8  m4 


and  R  =  -816. 50  x  10  "12  m6 

Since  the  coefficients  in  the  equation  are  real,  the  discriminant, 

3  2 

Q  +  R  ,  determines  the  nature  of  the  roots.  In  this  case, 

-24  1 2 

D  =  -726.2  x  10  m  ; 

so  all  the  roots  are  real  and  unequal.  S  and  T  are  then  computed  from 
Equations  (6.43)  and  (6.44)  and  are  displayed  in  Table  6.  From  Equations 

n 

(6.45),  (6.46)  and  (6.47),  the  solutions  for  u  *  a  are 

2 

u  j  =  1 . 309  x  10  4  m2 
u 2  =  -26.91  3  x  10-4  m2 
u3  =  0.953  x  10-4  m2 

Table  5.  The  Key  Variables  In  The  New  Theory 


Svmbo 1 
K 
L 
M 
N 
V 
P 


Value 

0.2669 

m  r 

/  S 

0.5435 

X 

,  n5  -1 
1 0  s 

0.2264 

X 

104  m/s 

24.65 

X 

in-4  2 
10  m 

-59.62 

X 

,  n-8  4 

10  m 

33.56 

X 

, n-l 2  ■ 

10  m 

Equation 

(6.42) 

(6.43) 

(6.44) 

(6.47) 

(6.48) 
(6.50) 


’*  V "V V  V  v  V  V  V  V  v\  V  V.V.V.V  A’A*aV\,Vv 
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Table  6.  Solutions  To  The  Cubic  Equation  Defining  The  Mass  Partition  In 

Conical  Flow 


9.348 

X 

10-4 

2 

m 

(cos 

59.37° 

9.348 

X 

IQ"4 

2 

m 

(cos 

179.37 

9.348 

X 

IQ'4 

2 

ra 

(cos 

299.37 

9.348 

X 

10"4 

2 

in 

(cos 

300.63 

9.348 

X 

10-4 

2 

m 

(cos 

180.63 

9.348 

X 

10"4 

2 

m 

(cos 

60.63° 

Since  u  must  satisfy 


o  <  u  <  1.046  x  10  4  ra2  , 


=  0.953  x  lO"4  m2 


The  radius  of  the  slug  is  9.76  mm  and  the  mass  division  is  such  that 

*2 


a(rC  +  rE) 


=  0.9101 


The  average  jet  velocity  is  the  sura  of  the  stagnation  point  velocity  of  the 
central  streamline,  3.598  km/s,  and  the  average  flow  velocity  into  the  jet 
in  the  stagnation  point  reference  frame,  obtained  from  Equation  (6.9),  3.232 
km/s.  The  result  is  an  average  jet  velocity  of  6.830  km/s.  The  average 
slug  velocity  is  the  difference  between  the  reference  frame  velocity  and  the 
average  flow  velocity  into  the  slug  obtained  from  Equation  (6.10),  2.969 
km/s.  The  result  is  0.629  km/s. 

Note  that,  in  this  case,  the  incoming  kinetic  energy,  in  the  stagnation 
point  reference  frame,  is  131.8  kJ/cra. 


IX.  COMPARISON  WITH  EXPERIMENT 

The  small,  systematic  deviations  between  theory  and  experiment 
discovered  by  Eichelberger*  will  now  be  discussed.  His  experimental 
observations  were  of  liner  mass,  jet  velocity,  and  slug  mass.  The  liner  -as 


♦Reference  4,  p.  402. 
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divided  into  sections,  and  the  slugs  corresponding  to  different  liner 
segments  were  recovered  by  firing  the  shaped  charge  into  a  reservoir  of 

water.  The  results  of  such  experiments  with  a  particular  charge  are  shown 

in  Figure  11.  The  jet  mass  was  estimated  by  subtracting  the  observed  slug 

mass  from  the  observed  cone  segment  mass. 

The  collapse  angle  as  a  function  of  liner  segment  position  was 
estimated  by  means  of  the  slug  mass  measurements  and  the  relation 


cos'  (?/2) 


(9.1) 


where  m  is  the  mass  of  the 
s 

segment.  Equation  (9.1)  is 

experimental  curve  of  m  vei 
r  s 

to  a  polynomial  which  was  si 

of  /J  displayed  in  Figure  12. 

the  cone,  along  the  axis  of 

erected  to  the  liner  wall. 

of  the  slug  is  no  longer  det 

collapse  angle. 


slug  and  m  is  the  mass  of  the  parent  liner 
the  differential  form  of  Equation  (3.5).  The 

8 

rus  m  (Figure  11)  was  fit  using  Birge's  method 

ubsequently  differentiated  to  obtain  the  values 
.  The  coordinate  x  is  measured  from  the  apex  of 
symmetry,  to  a  point  at  which  a  perpendicular  is 
Note  that,  according  to  the  new  theory,  the  mass 
termined  solely  by  the  mass  of  the  liner  and  the 


The  rapid  variation  in  /3,  especially  near  the  base  of  the  liner, 
indicates  that  the  collapse  of  this  shaped-charge  is  non-steady.  Pugh, 

Eichelberger  and  Rostoker*,  PER,  deduced  a  formula  for  the  shape  of  a 

conical  liner  collapsing  with  a  velocity  gradient  along  its  length  (instead 
of  through  its  thickness).  The  determination  of  the  jet-formation 
conditions,  the  collapse  velocity  V  and  the  collapse  velocity  gradient, 
dV  /dx,  was  relativaly  complicated,  but  the  jet  formation  itself  was  treated 
by  the  steady-state  theory.  Therefore,  any  discrepancies  found  between  the 
predictions  of  the  non-steady  jet-formation  theory  and  experiment  may  be  due 
to  either  inaccuracies  in  the  determination  of  the  initial  conditions  for 
jet  formation  or  inaccuracies  in  the  treatment  of  the  jet-formation  process 
itself.  The  collapse  velocities  and  jet  velocities  deduced  by  Eichelberger 
are  shown  in  Figure  13. 

The  nature  of  the  discrepancies  between  theory  and  experiment 
discovered  by  Eichelberger  is  illustrated  in  Figures  14  and  15.  In  Figure 
14,  the  theoretical  cumulative  mass  distribution  in  the  jet  (as  a  function 
of  velocity)  is  compared  to  the  results  of  several  experimental  mass 
determinations.  In  these  tests,  the  jet  first  perforated  a  target  plate 
then  its  tip  was  photographed  with  a  high-speed  camera,  and  finally  the 
remaining  jet  was  collected  and  the  particles  weighed.  The  jet  tip  velocity 
was  determined  from  the  photographs.  The  curve  was  computed  from  slug  mass 
data  using  the  "convergenceless"  PER  theory,  assuming  instantaneous 
acceleration  of  the  liner  elements  by  the  detonating  explosive  and 
deflection  of  the  liner  segment  according  to  the  Taylor  angle  hypothesis. 


*Reference  2,  p.  53b 
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Figure  13.  Eichelberger ' s  computed  values  of  collapse  velocity  V 

o 

velocity  V  as  functions  of  initial  position  of  liner  element 

(from  Reference  6) . 


In  Figure  15,  the  theoretical  spatial  distribution  of  the  jet  veloci¬ 
ties  is  displayed  along  with  the  experimental  observations  at  a  short  time 
( 3  5/x  s )  after  the  detonation  waves  passed  over  the  apex  of  the  liner.  The 
position  of  the  jet  tip  relative  to  the  detonation  wave  was  used  to  fix  the 
position  of  the  tip  for  the  theoretical  curve.  The  experimental  curve 
"leads"  the  theoretical  curve,  and  the  disagreement  is  greatest  at  the  rear 
of  the  jet. 

The  jet-mass  discrepancy  may  be  due  to  the  neglect  of  the  convergence 
gradient  in  the  treatment  of  the  steady-state  jet  formation  process.  A 
qualitative  explanation  of  the  difference  between  the  PER  and  the  observed 
jet  mass  distribution  may  be  obtained  from  the  new  theory.  The  differences 
between  the  jet  properties  calculated  by  the  Birkhoff,  MacDougall,  Pugh  and 
Taylor  (BMPT)  theory,  using  the  collapse  velocity  and  collapse  angle  of  the 
central  streamline  of  the  liner,  and  the  jet  properties  calculated  by  the 
new  theory  taking  into  account  a  linear  convergence-velocity  gradient  of 

strength  0.1107  km  s  '  ran  *  are  illustrated  in  Table  7.  For  these  "appar¬ 
ently  identical"  initial  conditions,  the  new  theory  predicts  a  jet  element 
with  about  8  percent  less  mass  and  about  3  percent  higher  velocity  than  that 
of  a  jet  element  whose  properties  were  estimated  by  the  BMPT  theory  which 
ignores  the  convergence-velocity  gradient.  For  the  slug,  the  new  theory 
predicts  an  element  slightly  more  massive  (less  than  1  percent)  and  13 
percent  faster  than  that  estimated  by  the  BMPT  theory.  If  a  correction  of 
92  percent  in  mass  and  103  percent  in  velocity  is  applied  to  the  curve  of 
Figure  19,  the  new  curve  is  in  better  agreement  with  the  experimental  data. 
This  is  shown  in  Figure  16.  Of  course,  there  is  no  reason  to  expect  that 
the  correction  would  be  constant  over  the  collapse  of  the  entire  liner,  but 
it  would  require  a  calculation  of  the  complete  jet-formation  process  to  find 
the  actual  corrections. 


An  explanation  of  the  difference  between  the  PER  theoretical  predic¬ 
tions  and  the  observed  jet  position,  probably  requires  more  than  just  the 
new  jet  formation  theory.  In  Figure  17,  the  point  J  on  the  axis  of  symmetry 
of  the  conical  liner  at  which  the  jet  is  formed  from  a  liner  point  origi¬ 
nally  at  P  is  controlled  by  a,  the  half-angle  of  the  cone  apex  and  by  the 
deflection  angle  5.  For  steady-state  collapse,  the  angle  8  is  the  Taylor 
angle  used  in  the  PER  theory  for  the  case  of  the  unsteady  collapse.  For 
liner  collapse  in  which  there  is  a  velocitv  gradient  along  the  length  of  the 
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liner,  as  in  the  case  of  unsteady  collapse,  Chou,  Hirsch  and  Ciccarelli  , 

CHC ,  have  shown  that  the  deflection  angle  is  no  longer  the  Taylor  angle,  but 
depends  on  the  collapse  velocity  gradient.  Thus,  both  the  collapse  angle  ft 
and  the  deflection  angle  8  depend  on  the  collapse  velocity  gradient  along  the 
length  of  the  liner. 


According  to  some  calculations  made  by  CHC  for  the  collapse  of  the 
BRL  81mm  shaped  charge,  the  difference  between  the  deflection  angle 
"observed"  in  computer  simulations  of  the  cone  collapse  and  the  Taylor  angl 
increases  from  some  13  percent  near  the  apex  to  as  much  as  75-80  percent 
near  the  base  of  the  cone.*  Similar  results  would  apply  to  the  charge 
analyzed  by  El chelberger  ,  so  it  appears  that  the  deflection  angle  was 
probably  underestimated  by  PER  and  so  the  distance  of  the  collapse  point  J 
from  the  cone  apex  was  probably  underestimated.  The  combination  of  revised 
input  conditions  and  the  slightly  higher  jet  velocities  obtained  from  the 
new  theory  may  account  for  the  position  discrepancy  at  short  times  after 
forma  t ion . 


Table  7.  COMPARISON  OF  THEORETICAL  RESULTS 


QUANTITY  IN  STAGNATION-POINT  FRAME 

THEORY-EQUATION 

1  # 

Incoming 

Flow  Velocity 

7  - 

vT  (9,15)  =  3.038  km/ s 

F 

BMPT- (3 . 1 ) 

= 

<Vrt,>  =  2.99  7  km/s 

FN 

NEW- (6. 8) 

I ncoming 

Flow  Kinetic  Energy 

KE.  = 
l  n 

135.4  kJ/cm 

BMPT- f  B . 3 ) 

= 

131.8  kJ/cm 

NEW- ( C . 3 ) 

3  . 

Outgoing 

Jet  Flow  Velocity 

v7i  = 

Vp  =  3.038  km/s 

BMPT- (3.  3') 

= 

'<V^>  '  =  3.232  km/s 

NEW- (6 .9) 

U  . 

Outgoing 

Jet  Kinetic  Energy 

KElout  - 

13.2  kJ/cm 

BMPT 

= 

14.2  kJ/cm 

NEW 

5. 

Outgoing 

Slug  Flow  Velocity 

i  v‘~k 

F2 

' V  1  =  3.038  km/s 

F 

BMPT-C  3.4) 

= 

=  2.96  9  km/s 

F2 

NEW- (6 .10) 

♦Reference  9,  p.  23. 


6.  Outgoing  Slug  Kinetic  Energy 


KE0  „  =  122.2  kJ/cm 

Zou  t 


BMPT 


=  117.4  kJ/cm 


NEW 


QUANTITY  IN  LAB  FRAME  FOR  Vgp  =  3.598  km/s 

7 .  Jet  Velocity 

V.  =  6.636  km/s 

J 

=  6.830  km/s 

8 .  Jet  Mass 

nr  =  0.0976  m^  =  2.865  g/cm 

=  0.0899  tn p  =  2.639  g/cm 

°.  Slug  Velocity 

V  =  0 . 5  6  0  km / s 

s 

=  0.629  km/s 

10.  Slug  Mass 


THEORY-EQUATION 

BMPT-( 3.3) 
NEW 

BMPT 

NEW 

BMPT-( 3.4) 
NEW 


n 


s 


0.9024  m^ 
0.9101  ffll 


26.491  g/cm 
26.717  g/cm 


BMPT 

NEW 


X.  SUMMARY  AND  CONCLUSIONS 


This  study  appears  to  be  the  first  attempt  to  explicitly  integrate 
convergence  effects  into  the  jet-formation  equations.  For  initial 
conditions  appropriate  to  a  conventional  conical-lined  shaped  charge,  the 
predicted  changes  to  the  results  of  the  previous  theory  are  small,  but  in 
the  proper  direction  to  explain  some  of  the  discrepancies  observed  by 
E i che 1  be rge r  . 


The  new  theory,  in  both  wedge  and  conical  flow  form,  conserves  mass, 
kinetic  energy  and  linear  momentum.  Appendix  A  demonstrates  the  conserva¬ 
tion  of  kinetic  energy  in  wedge  flow;  Appendix  B  does  the  same  for  "uniform 
conical  flow."  Appendix  C  illustrates  the  key  results  in  the  demonstration 
of  conservation  of  kinetic  energy  in  nonuniform  conical  flow  with  a  linear 
velocity  gradient.  For  the  rather  artificial  condition  of  constant 
energy  influx  and  varying  strength  of  the  convergence-velocity  gradient,  the 

mass  partitioning  for  the  wedge  flow  is  confined  between  mg/m  =  cos  2  (0/2) 

and  m  /m  =  cos  (3/2),  for  positive  values  of  the  minimum  collapse 

velocity  V  in  the  case  considered.  Very  similar  results  are  found  for  the  ca 

of  conical  flow. considered.  It  appears  that  the  strengths  of  the  convergence- 
velocity  gradients  for  conventional  conical  liners  lie  in  the  range  0.1  to 
0.°  !cn  s-’1  tin'-  . 

A  HEMP  simulation  of  the  collapse  of  the  shaped  charge  used  in 
Eichelberger ' s  experiments  and  a  calculation  of  the  integrated  effect  of  the 
convergence-velocity  perturbation  ;ould  be  required  to  show  whether  the 
convergence-velocity  effect  alone  is  sufficient  to  bring  theory  and 
experiment  into  agreement.  The  HEMP  code  simulates  the  collapse  dynamics  of 

a  shaped-charge  liner  rather  accurately  ^ J ,  but  cannot  properly  form  the  jet 
and  slug  for  most  shaped- charge  designs.  The  HEMP  input  would  be  used  in 
the  new  jet  formation  equations  since  it  "automatically"  includes  the 
unsteady  collapse  effects. 
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APPENDIX  A 


Conservation  of  Kinetic  Energy 
in  Wedge  ^lov 


KE  . 
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-w 

•'  o 


a  V2(x) 


dx 


-hf 

J  o 


(V  +  clx)  1  dx 
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o  a  O')  CX  q 

-a  +e—  a-  +  p-g  a 


out 
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V"  aV  2  ■? 
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where  a  +  a  =  a 


APPENDIX  B 


Conservat ion  of  Kinetic  Energy  in 

"Vniform  Conical  Flow" 

1  fa  V“ 

=  y  .:A  ~s.  2”  (r  -  cx)dx 

in  2  o  )  A  n 

J  O  o  w 

(B-l) 
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(B-4) 
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=  ~  ':o  a  fr_  "  ta/2! 
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APPENDIX  C 


Conservation  of  Kinetic  Energy  in  Nommifora 
Conical  Flow  With  Linear  Velocity  Gradient 


:  -i,»  f 

in  2  o  3 


a  V|N  (x) 


2^  ( r  -  ex)  dx 

C 


CVCFN  +  trC  '  '?X)dX 


f r CVCFNa  +  (2‘rC  C''  DFN5  VOFN  2 


(r^:  -  2caV0FJJ)  f-  -  au/4  ' 


The  fraction  of  the  area  of  the  incoming  flow  through  which  the  jet  material 
passes  is  illustrated  in  Figure  C-l.  The  area  intercepted  by  the  incoming 
jet  flow  is 


A  =  -ry'  (r  *  +  r  ,) 
y'  *  a*  V 


r  =  r 


=  r  -  c a’ 
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(2r  -  2f.a*  -  ry’) 
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going  get  flow,  v:  . 
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(C-9) 


and 
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(C-10) 
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(C-12) 


gives  the  streamline  spacing  relation.  Note  that,  when  y'  =  0,  x"  =  0,  and 


when  y'  =  a  -  a*,  X"  - 


[(r  -  za*)~  -  (r  -  ea)2]'172 

C  2  C 
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KE 


out 


i  [v;N 


2^x"dx" 


*  j  oA,  3 
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where  the  first  term  represents  the  energy  per  unit  lenqth  in  the  jet  and  the 
second  the  energy  per  unit  length  in  the  slug. 
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For  the  slug 
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(C— 14) 
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2  c  2  C  2 


(C— 16) 


For  the  jet. 
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Note  that  V'  (x')  decreases  as  x'  increases,  while  V"  (x  )  increases  as  x 
FN  FN 

increases. 


Substituting  Eqs.  (C-14)  and  (C-18)  into  Eq.  (C-13)  and  integrating  gives 
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AMXSY-GI  (B.  Simmons) 

Cdr ,  USATEC0M 

ATTN:  AMSTE-T0-F 

Cdr,  CRDC ,  AMCC0M 
ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-SPS-IL 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  comments/answers  to  the  items/questions  below  will 
aid  us  in  our  efforts. 


1  .  BRL  Report  Number _ 

2.  Date  Report  Received_ 


Date  of  Report 


3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or 
other  area  of  interest  for  which  the  report  will  be  used.) _ 


4.  How  specifically,  is  the  report  being  used?  (Information  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
as  man-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved, 
etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 


CURRENT 

ADDRESS 


Organization 


Address 


City,  State,  Zip 


7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the 
New  or  Correct  Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


OLD 

ADDRESS 


Organization 


Address 


City,  State,  Zip 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 
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